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ABSTRACT 

We report the first two-dimensional mapping of the stellar population and 
non-stellar continua within the inner 180 pc (radius) of NGC 1068 at a spatial 
resolution of 8pc, using integral field spectroscopy in the near-infrared. We have 
applied the technique of spectral synthesis to data obtained with the instrument 
NIFS and the adaptive optics module ALTAIR at the Gemini North Telescope. 
Two episodes of recent star formation are found to dominate the stellar popu- 
lation contribution: the first occurred 300 Myr ago, extending over most of the 
nuclear region; the second occurred just 30 Myr ago, in a ring-like structure at 
~ 100 pc from the nucleus, where it is coincident with an expanding ring of H 2 
emission. Inside the ring, where a decrease in the stellar velocity dispersion is ob- 
served, the stellar population is dominated by the 300 Myr age component. In the 
inner 35 pc, the oldest age component (age>2Gyr) dominates the mass, while 
the flux is dominated by black-body components with temperatures in the range 
700 < T < 800 K which we attribute to the dusty torus. We also find some con- 
tribution from black-body and power-law components beyond the nucleus which 
we attribute to dust emission and scattered light. 

Subject headings: galaxies: individual(NGC 1068) — galaxies: active — galaxies: 
nuclei - galaxies: stellar population 
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Introduction 



NGC 1068 is the closest, brightes t and for this reason, the m ost studied Seyfert 2 



galaxy. Its distance of only 14.4 Mpc ( ICrenshaw &: Kraemer 



20001 ) corresponds to a scale 



of 70 pc arcsec 1 at the galaxy. In spite of its proximity, there is only limited information 
about the stellar population in the vicinity of its active nucleus. 



Lynds et al.l (119911 ) reported a high contrast structure surrounding the nucleus in a 
Hubble Space Telescope (HST) Wide Field and Planetary Camera 2 (WFPC2) optical 
image of NGC 1068 in the continuum (centered at 5470 A) interpreted as a stellar cluster. 



Thatte et al 



( 11997a! ) reported the detection of this stellar cluster from stellar CO absorption 



features in the H and K b and spectra and estim a ted it s size as FWHM« 50 pc and age in 



the range 5 — 16 x 10 8 yr. 



Crenshaw fc Kraemer! (120001 ) have analysed two long-slit spectra 



obtained with the HST Imaging Spectrograph (STIS) covering the inner ~500pc in the 
wavelength range 1150-10270 A. They have found that the continuum in this inner region 
has contributions from two components: electron-scattered light from the hidden nucleus 
and stellar light. They conclude that the electron-scattered light dominates in the UV, 
while stellar light dominates at optical wavelengths. They found that the scattered light is 
strongest in a hot spot at ~ 0'.'3 northeast of the nucleus, being enhanced also in regions 
of strong narrow-line emission. The stellar component was found to be more concentrated 
southwest of the hot spot, approximately coincident with the nucleus, being dominated by 
an old (age > 10 9 yr) stellar population. 



Crenshaw fc Kraemer! (120001 ) also support the presence of a nuclear stellar cluster 
~ 200 pc in extent and argue that it cannot be too young (in the sense of having blue 
ionizing stars), because the light conce ntration seen in the optical continuum is not 



seen in the UV (IMacchetto et al. 



Macchetto et al. 



1994J). The UV Faint Object Camera HST images of 



19941 ) resemble instead the narrow-band [OIII] optical images, consistent 
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with the interpretation that they are due to scattering, which is also supported by the 
fact that this light is highly polarized. As is well known by the AGN community, it 
was the discovery of this polarized light emission, showing; a Seyfert 1 spectr um, that 



led to the Unified Model of AGN f Antonucci fc Miller 



1985 



Antonucci 



19931 ). Another 



argument against the presen ce of ionizing stars is the absence, in the UV spectra of 



Crenshaw fc Kraemerl (120001 ). of stellar absorption lines characteristic of the atmosphere of 



such stars. 



More recently, lEmsellem et al.l (120061 ) . using optical observations with the SAURON 
integral field spectrograph, have claimed that this galaxy is a candidate for a stellar velocity 
dispersion (er*) "drop" in the central ~2 arcseconds: a* decreases to ~100kms _1 in this 
region, being surrounded outwards by higher values in the range 150-200 km s -1 . This 
<T*-drop has been interpreted as the result of central gas accretion followed by an episode of 
star formation. 



Davies et al.l ( 120071 ) have also observed NGC 1068 and other nearby AGN using the 
near-IR integral field spectrograph SINFONI at the Very Large Telescope (VLT). In most 
of the galaxies they found signatures of circumnuclear starbursts with ages of 10-300 Myr 
(i.e., no longer ionizing). In the case of NGC 1068, they confirm the presence of a ovdrop, 
and report also an excess of stellar continuum within the inner 1", which they attribute to 
the presence of a nuclear disk that has experienced recent star formation. This would also 
explain the decrease in velocity dispersion. 



Martins et al 



(120 lOf ) performed stellar population synthesis on the nuclear and 
extended regions of NGC 1068 by means of near-infrared spectroscopy (0.8 — 2.4/im) 
using the NASA Infrared Telescope Facility SpeX-XD data. A 0'/8 x 15" slit oriented 
in the north-south direction was employed. Their main result is that traces of young 
stellar population are found at ~100pc south of the nucleus. They also report a ~25% 
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contribution of a power-law continuum at the nucleus and surrounding regions, with excess 
contribution (up to 60%) of this component at 100-150 pc to both sides of the nucleus, 
associated by them with regions where new stars are being formed. They also report that 
hot dust (1000K<T< 1400 K) has an important contribution to the nuclear region while 
cold dust appears mostly to the south. A significant contribution of an intermediate-age 
stellar population to the continuum, especially in the inner 200 pc, was also found. 

With the goal of mapping and quantifying the processes of feeding and feedback 
in the active nucleus of NGC 1068, as well as its stellar population, we have analyzed 
two-dimensional spectroscopic data obtained with the Gemini North Near- Infrared Integral 
Field Spectrograph (NIFS) in the J, H and K bands. In the present paper, we report the 
mapping and analysis of the stellar population and near-infrared continuum within the 
inner 180 pc radius of NGC 1068 at a spatial resolution of ps 8 pc. We also include a brief 
analysis of the stellar kinematics, with the goal of relating it to the stellar population 
characteristics. 

This paper is organized as follows: in Section 2 we discuss the observations and data 
reduction procedures, in Section 3 we derive the stellar kinematics, in Section 4 we obtain 
the properties the stellar population, in Section 5 we discuss the results, comparing the 
stellar population properties with the derived kinematics and the molecular gas distribution 
and in Section 6 we present our conclusions. 

2. Observations, data reduction and analysis 

Spectra of NGC 1068 were obtained with the Near-infrared Integral Field Spectrograph 
(NIFS; McGregor et al. 2003) on the Gemini North telescope in December 2006 under 
program GN-2006B-C-9. NIFS was used with the ALTAIR facility adaptive-optics system 
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in its natural guide star mode. The compact Seyfert nucleus of the galaxy was used as the 
adaptive-optics reference object. A standard star observed under the same conditions as 
NGC 1068 has a spatial profile with a FWHM of ft' 12 at both H and K, corresponding to a 
spatial resolution of ~ 8 pc at the galaxy. 

NIFS is an image-slicer integral-field spectrograph with a square field of view of 3"x3", 
divided into 29 slitlets each 0'.'103 wide with a sampling of ft'044 along each slitlet. The 
instrument was set to a position angle on the sky of 300° to align the slitlets approximately 
with the axis of the radio jet and ionization cone. This provided coarser spatial sampling 
along the jet axis, and finer spatial sampling across the jet. 

Observations at the Ks grating setting were obtained on 12 December 2006 (UT). 
Similar observations at the Kl grating setting were obtained on 13 December 2006 (UT), 
in the H band on 14 December 2006 (UT), and in the J band on 15 December 2006 (UT). 
The same orientation, observing procedure, exposure time, and standard stars were used on 
each of these nights. The spectra used in the present study cover the H band from 1.48 /im 
to 1.79 /im and the Kl band from 2.1 /im to 2.5 /im, both at a spectral resolving power of 
~ 5290. The spectra in other wavelength bands will be discussed in future papers. 

Each data set consisted of multiple 90 s exposures. An offset sky position was recorded 
first and was followed by nine galaxy fields obtained on a 3 x 3 frame grid centered on 
the galaxy nucleus. The frame offset was l'.'O, so frames were recorded at offsets of — l'.'O, 
ft.'O, and + l'.'O with respect to the nucleus along and perpendicular to the slitlets. This 
resulted in a full field of view (FOV) 5'.'0 x 5'.'0 on the sky with maximum exposure in 
the central 3'.'0 x 3'.'0 and less exposure in the corners of the full field. The exposure grid 
was repeated nine times at each grating setting for a total exposure time on any one of 
the nine fields of 810 s. No pixel-level dithering of the field positions was employed. The 
nucleus of NGC 1068 was within each field, and so was subsequently used to spatially 
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register each exposure. An argon/xenon arc exposure was obtained immediately after the 
galaxy exposures were completed. The A2 telluric standard star HIP 5886 was observed 
immediately prior to the NGC 1068 observation, and the AO star HIP 18863 was observed 
immediately after the NGC 1068 observation. Both were selected so that they were observed 
through approximately the same atmospheric column as the NGC 1068 observations. 

The full FOV and orientation of our NIFS observations are illustrated by the white 
square over an Hubble Space Telescope (HST) image through the filter F606W in the 
left panel of Fig.fU where we also show, for reference, within our FOV, a narrow-band 
[OIH1A5007A image (top right panel of Fig. [I]) showing the known ionization cone 



( iSchmitt et al. 



20031 ). as well as the radio jet (bottom right panel). 



The data reduction was performed usin g the GEMINI IRAF N 



followed standard procedures, as described in 



Storchi-Bergmann et al 



FS pa ckage and 



fl2009[ ). The nine sky 



exposures were first median-combined and dark subtracted to form an average sky frame. 
Individual galaxy exposures were then dark subtracted and the average sky exposure was 
subtracted. The 29 two-dimensional slitlet spectra for each exposure were cut from the 
raw data image and flatfielded. Known bad pixels identified from the flatfield and dark 
exposures were corrected by linear interpolation. Each two-dimensional slitlet spectrum was 
then transformed to a rectilinear coordinate grid with horizontal wavelength and vertical 
spatial axes using transformations derived from the arc lamp and spatial calibration mask 
exposures. These transformed spectra were then stacked in the second spatial dimension to 
form a three-dimensional data cube with two spatial and one spectral axes. 

The individual data cubes were then registered spatially by collapsing each data 
cube in the spectral direction to form a continuum image, determining the centroid of 
the NGC 1068 nucleus in each image, and then spatially shifting and median-combining 
all galaxy cubes using the IRAF imcombine task. The resulting mosaiced data cube 
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Fig. 1. — The large left panel shows an HST image of NGC 1068 through the filter F606W 
with the FOV of the NIFS observations indicated by the square. The top righ t panel shows 



Schmitt et al. 



(120031 ) . as observed 



an intensity map of the [OIII]A5007 emission line from 
within the same FOV of the NIFS observations. The bottom right panel shows the radio jet. 
The black cross marks the position of the galaxy nucleus. 
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was corrected for telluric absorption by multiplying each spectrum within the cube by a 
one-dimensional correction spectrum derived from the one-dimensional spectrum of the 
telluric standard star, HIP 5886, extracted from its data cube. A similar large-aperture 
one-dimensional spectrum of the total light of HIP 5886 was then used to flux calibrate the 
mosaiced galaxy data cube. The intrinsic spectrum of the standard star was set to the level 
of its 2MASS magnitude in the wavelength band of the spectrum, and its shape was set to 
that of a 8600K blackbody, which has the same 2MASS J - K color as HIP 5886. 

In order to have approximately square spaxels, we have then binned together two pixels 
together along the slices. The resulting calibrated data cube contains 2376 spectra covering 
an angular region 4'/5 x 5'/0 (along the x axis and y axis, respectively) at a sampling of 
0'.'103 x 0'.'092, which corresponds to ~ 315 x 350 pc 2 (or a radius of ~180pc about the 
nucleus) and ~ 7.0 x 6.4 pc 2 , respectively, at the galaxy. 



2.1. Data analysis 

2.1.1. Line-of-Sight Velocity Distributions 

In order to map the stellar kinematics, we have obtained the line-of-sight velocity 
distributions (LOSVD) by fitting the absor ption spectra in the H and K bands using the 



penalized Pixel-Fitting (pPXF) me t hod o f 



procedure described in 



Riffel et al 



Cappellari fc Emselleml (120041 ). following the 



( 120081 ). In summary, the pPXF method finds the best 
fit to a galaxy spectrum by convolving st ellar template spectra with a given LOSVD | X(t;)|, 
repre sented by Gauss-Hermite series (e.g. 
2010): 



van der Marel & Franx 


1993: 


Gerhard 


1993; 


Riffel 



L(v) 



-(V2)/y 2 



a*v27T 



M 



1 + 5^ h m *H m (y) 



m=3 



X 



- 10 - 



where y = (v — V*)/a, V* is the stellar centroid radial velocity, a* is the velocity 



dispersion, v = c In A and c is the spe ed of light. The H m are the 



h m * are the Gauss-Hermite moments (jCappellari fc Emsellem 



iermite polynomials, 



20041 ) . The pPXF routine 



outputs the stellar centroid velocity (K), velocity dispersion (a*) and t h e hig her 



order Gauss-Hermite moments 



and h i*. As discussed in 



Riffel et al 



(120081 ) and 



Winge. Riffel fc Storchi-Bergmannl (120091 ). the use of a library of stellar templates, instead 
of a single stellar spectrum, is fundamental for a reliable derivation of the stellar kinematics. 
In the K-band, we have used as templates the spectra of 60 stars from the the Gemini 
library of late spectral ty pe stars observed with the Gemini Near- Infrared Spectrograph 



(GNIRS) IFU and NIFS (IWinge. Riffel k Storchi-Bergmannl l2009h . In the H-band, we do 



not have a large number of templates available, and we have thus used the spectra of five 
late type stars with public NIFS data in the Gemini archive. The H-band spectra of these 
stars are shown in Fig.[2l together with the NGC 1068 nuclear spectrum. 



2.1.2. Spectral Synthesis 



In order to map the distribution of the continu um and ste 



have u sed the starlight spectral synthesis code (lAsari et a 



2005b 



aj) adapted to spectra in the near-IR (IRiffel et al 



ar population properties, w e 



2007 



Cid Fernandes et al. 



2009). This code fits a galaxy's 



spectrum by a combination of elements from a "spectral base" (desc ribed below), giving a s 



output the weights of e ach base element in the synthetic spectrum ( iCid Fernandes et al. 



2009 



Riffel et al. 



20091 ). 



The spectral ba se was constructed with the m ost recen t Evolutio n ary Populat ion 



Synthesis models of 



following 



Riffel et al 



Maraston fc Stromback! ( 2011 



see also 



Maraston! ( 1998 



20051 )) 



(120091 ). and comprise single stellar populations (SSPs) synthetic 



spectra covering 12 ages (t =0.01, 0.03, 0.05, 0.1, 0.3, 0.5, 0.7, 1, 2, 5, 9 and 13 Gyr) 
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HIP1 13543 G2V 



HD35833 GO 



HD40280 KOI 1 1 



HIP1 1 1925 KOI 



HIP601 45 MO 







NGC1068 Spectrum 








1.50 1.55 



1.60 1.65 
A (/im) 



1.70 1.75 



Fig. 2. — H-band spectra of the stars used as templates by the pPXF, as well as sample 
spectrum of NGC 1068 extracted at 0'.'5 NE from the nucleus (position A from Fig.[3]). 
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and solar metallicity. We did not include a stellar population younger than lOMyr as we 
verified it was not necessary b ecause it becomes degenerated with a power-law component 



( jStorchi-Bergmann et al 



20001 ) also included in the synthesis. The power-law (F v oc v 



was used to repres ent the possible contributi on of a featureless continuum from the 



active nucleus (e.g. 



Cid Fernandes et al. 



2004J ) . We also included black-body functions for 



temperatures in the range 700-1400 K (one for each 100 K interval) in order to account for 
possible contributions from dust emission. 



Following 



Cid Fernandes et al 



( 120041 ). we have binned the contribution of the SSP Xj 
elements into stellar population components of four age ranges tf young (x y : t y < 50 Myr); 
young- intermediate [x y i. 0.1 < t yi < 0.7 Gyr), intermediate-old (xi Q : 1 < t io < 2 Gyr) and 
old (x Q : 5 < t Q < 13 Gyr). We present the synthesis results in two forms: percent flux 
contribution of each SPC at A = 2.1/xm (xj) and percent mass contribution to the total 
mass (mA 



Results: Stellar Kinematics 



As de scribed in the previous section , we have used the penalized Pixel-Fitting (pPXF) 



method of 



Cappellari fc Emselleml (120041 ) to fit the absorption spectra in the H and K bands 



in order to obtain the 
(e.g. 



Riffel et al. 



2008 



ine-of-sight velocity distributions (LOSVD). In our previous studies 



2010 



20111 ) we have used only the K-band spectra for the derivation 



of the stellar kinematics. But in the case of NGC 1068, the K-band spectra within w 0'.'3 
from the nucleus are strongly diluted by a red continuum that can be attributed to emis sion 



by hot dust from the circum-nuclear torus, as in the case of NGC 4151 ( jRiffel et al 



2009|). In 



addition, the stellar K-band ab sorption features (CO 



the ionization cone is located (j Crenshaw fc Kraemer 



jand s) towards the north-west, where 



20001 ). are contaminated by strong 



emission from the [Ca vm]A2.322 /zm coronal line. In order to avoid these contaminations 



13 



present in the K-band spectra, we have used the H-band spectra to derive the stellar 
kinematics for angular distances from the nucleus smaller than 0'.'3. For the other regions 
we used the K-band spectra whenever the contamination by the [Ca vm]A2.322 /mi coronal 
line was negligible, as the stellar absorption features are stronger in the K-band than in the 
H-band, resulting in a better constrained stellar population, particularly for the outermost 
regions. In order to "automatically" select which band to use for the derivation of the 
stellar kinematics, we used the flux of the [Ca vm]A2.322 fim coronal line: in regions with 
fluxes of this line higher than 3a (a here meaning the root mean square deviation of the 
continuum flux in the neighboring continuum), we used the H-band spectra, while in other 
regions we used the K-band spectra. 

In the top and bottom panels of Figure[3j we show representative fits (in red) compared 
to the observed spectra (in black), extracted at locations identified in the resulting velocity 
dispersion map: the nucleus, 0'/5 north-west, l'.'O south-east and l'/5 north-west of the 
nucleus. The region where we have used the H-band spectra for the derivation of the 
kinematics is shown in the center right panel of Fig.[3] (shaded region in the figure). 

The left panel of Fig. H] shows the centroid velocity map, from which we subtracted 
the systemic velocity of Vlsr = 1168 ± 15kms _1 [corrected for the motion relative to the 
local standard of rest (LSR)] adopted as the average stellar velocity within an aperture of 
0'/5 x 0'/46 centered on the nucleus. The centroid velocity map shows a rotation pattern, 
with the kinematic major axis approximately along the east-west direction, at PA^ 80°, 
with blueshifts to the east and redshifts to the west, reaching a projected velocity amplitude 
of 45kms~ 1 . One-dimensional cuts showing the centroid velocities along two "pseudo-slits" 
with width 0'/3 oriented along the galaxy major and minor axes (PA?^ 80° and ~ —10°, 
respectively) are shown in Fig.[5j The pseudo slit along the major axis is illustrated in the 
central left panel of Fig.[3]by the dashed lines. 
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Fig. 3. — Top and bottom: sample spectra showing the fits obtained via the Penalized Pixel 
Fitting method (in red) compared to the observed spectra (in black) at locations identified 
in the central panels. Center left: velocity-dispersion map showing the pseudo-slit along the 
major axis. Center right: "mask" showing the region where the kinematics were derived 
from the H-band spectra. 
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The right panel of Fig.H] shows the velocity dispersion a* map. The highest values 
reach o+ ~ 190 km s -1 , observed right at the nucleus and in patchy regions farther than 1" 
from the nucleus towards the borders of the field. Within a region between « (X'2 and ~ 1" 
from the nucleus (corresponding to 14 pc and 70 pc, respectively, although more extended 
along east-west than along north-south), a* decreases to 125 km s -1 . One-dimensional cuts 
showing the velocity dispersions along the two "pseudo-slits" - passing through the galaxy 
major and minor axes are shown in Fig.|5j 

In order to calculate the errors in the centroid velocity and velocity dispersions, we 
have performed Monte Carlo simulations in which Gaussian noise was added to the spectra. 
The corresponding error bars are illustrated in the one dimensional cuts of Fig.[5j showing 
that the errors both in the centroid velocities and in the velocity dispersions range from 10 
to 20 krns" 1 . 



As described in Sec.|2 



to spectra in the near-IR ( IRiffel et al 



4. Results: Stellar population 



we have used th e starlight spectral synthesis code adapted 



20091 ) to fit the galaxy's spectra by a combination 
of templates from a base including black-body and power-law continua and single age 
stellar pop ulations, giving as output the contribution o f each base element to the synthetic 



spectrum (jCid Fernandes et al. 



2009 



Riffel et al 



20091 ) 



We show in Fig.[6j the results of the synthesis for two spectra (extracted within our 
aperture of / /lx0 / . / 09), one at the nucleus (top panel) and the other at 1" north of it 
(bottom panel). The observed (continuous line) and synthetic (dashed line) spectra were 
normalized to unity at 2.1 /im. 



The spatial distribution of the stellar population components are shown in Fig. [71 
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Fig. 4. — Stellar kinematics obtained from the fit of the stellar population features in the H 
and K bands (see Fig.[3]). Left: centroid velocity map; right: velocity dispersion map. The 
central cross marks the position of the nucleus and the color bars show velocity values in 
kms" 1 . The orientation is indicated by arrows in the left frame. 
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Fig. 5. — Stellar kinematics extracted along two pseudo-slits shown together with typical 
error bars. Left: centroid velocity; right: velocity dispersion. Top: along the galaxy minor 
axis (PA= —10°); bottom: along the major axis (PA=80°). 
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1 .5 2 2.5 

A(ytim) 



Fig. 6. — Fit of the combined stellar population and continuum components to the nuclear 
spectrum (top) and to an extra-nuclear one at 1 "(70 pc) north of the nucleus (bottom). The 
observed spectra are shown as continuous lines and the fits as dashed lines. The nuclear 
spectrum is dominated by the contribution of the blackbody functions corresponding to the 
dust emission. 
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The top panels show the percent contributions to the light at A2.1yum of the young, 
young-intermediate, intermediate-old and old age components, from left to right: x y , x yi , 
x io and x Q , while the bottom panels show the corresponding mass contributions m y , m yi , 
m.i and m G . We have masked out from the figures the results of the synthesis in the regions 
for which adev (see below) is smaller than 2% in order to show only the results based on 
robust fits. 

The flux contribution to the continuum within the inner 180 pc of NGC 1068 is 
dominated by the young and intermediate age stellar populations. The flux contribution 
of the young population x y , as well as its contribution to the total mass m y , shows a 
"hole" of radius ~ 1" around the nucleus covering the east-north-west directions, but 
extending almost to ~2" towards the soutwest-south. The x y and m y are distributed 
around this "hole" , with enhanced contributions in a partial ring around the nucleus from 
the northeast-east towards the south and to the west-southwest. The young- intermediate 
age contribution x yi and m yi shows also a "hole" around the nucleus but with approximately 
half the radius of the one observed in x y and m y , and with similar extent in all directions. 
The southwest-south extension of the "hole" seen in x y and m y is not present in x yi and 
m y i, which show enhanced contributions in this region - a clump between rj 1 and 2" to 
the right of the nucleus in Fig. [7]- where m y i w 100%. The remaining regions of enhanced 
contributions of x yi and m yi seem to be anti-correlated with those of enhanced x y and m y , 
being observed mostly to the north-east through the north and north-west. 

The contribution of the old stellar population, x Q and m Q , is restricted to within ^0'.'5 
(35 pc) from the nucleus, where the contribution to the total mass m is dominated by this 
component. The contributions of the intermediate-old stellar population Xj G and are 
small and avoid the nuclear region, being somewhat co-spatial with the regions of enhanced 
contribution of the young stellar population. 
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Fig. [8] shows the contribution of the power-law and black-body components, together 
with the reddening derived from the stellar population synthesis. The rightmost panel 
shows the distribution of the adev, which is the percent mean deviation of the fit - 



Q,-M,\/Q,x 100 



(jCid Fernandes et al 



whe r e Q\ is the observed spectrum and M\ is the fitted model 



2004 . 



2005bl ) - and gives the robustness of the fit. Values of 
adev < 2% correspond to good fits, while higher values, observed towards the corners of 
the field, indicate less robust fits. 

Fig. |H] shows that the contribution to the flux within 0'.'5 of the nucleus is dominated 
by the black-body components with 700 < T < 800 K, which also show a significant 
contribution to a radius ~ l".The black-body components with 900 < T < 1000 K also have 
some contribution closer to the nucleus, while for T>1100K the contribution occurs even 
closer to the nucleus but is of negligible strength. Fig.E] also shows that the reddening 
obtained from the synthesis is only significant within the inner w 0'.'5, where it reaches 
E(B-V) = 1.5 mag. The power-law component is extended towards the south-southeast (right 
in Fig. El), while the vertical structure is an artifact probably due to scattering of the strong 
nuclear source along the NIFS mirror slices. 



5. Discussion 



5.1. Stellar kinematics 



Fig.H] shows a strong signature of rotation about the nucleus and Fig.0 reveals a 
steep rotation curve, with the line-of-sight velocity reaching ^40kms _1 at just 0?2, or 
140 pc from the nucleus. This result can be compared with those of previous studies. The 
first two-dim ensional study of the s tellar kinematics of the nuclear region of NGC 1068 



was done by 



Garcia-Lorenzo et al. 



(119971 ) using the instrument 2D-FIS on the William 
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Fig. 7. — Top: spatial distributions of the percentage contribution of each SPC to the flux 
at A = 2.1 /iin. From left to right: young (< 100 Myr), young-to-intermediate (0.3-0.7 Gyr), 
intermediate-to-old (1-2 Gyr) and old (5-13 Gyr). Bottom: spatial distributions of the 
mass-weighted contribution from each SPC. 
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Fig. 8. — Top, from left to right, distributions of: (1) percent flux contributions of the 
power-law component; (2)reddening E(B-V) derived from the stellar population synthesis; 
(3) the parameter Adev, which quantifies the quality of the spectral fit. Bottom, from left to 
right: contributions of the blackbody components of temperature ranges shown at the top 
of each panel. 
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Herschell Telescope (WHT). 95 fibers, each of 079 diameter, were used to cover the inner 
972 x 1272 at 5A spectral resolution in order to map the stellar kinematics in the Call 
triplet around 8500A. The main result of this paper was that the stellar kinematics of the 
inner 3" has a center which does not agree with the center of the stellar kinematics at 
radii larger than 5", suggesting a decoupled kinematics between the the inner and outer 
regions. Our observations are restricted only to the inner region, thus we cannot verify the 
claimed decoupling. But in the inner region, the velocity amplitude of ps50kms -1 , and the 
orientation of the major axis, of 88° ± 5°, are in agreement with our velocity field shown in 
Fig.|H Their velocity dispersions - in the range 250-300 km s _1 - are somewhat larger than 
our values (120-180kms _1 ). This difference is probably due to the different methods and 
stellar template spectra used to obtain the velocity dispersion values. 



We can also compare our results with those from iDavies et al.l (120071 ) , who observed the 
inner « 2" of NGC 1068 at 071 resolution in the H-band with SINFONI at the VLT. They 
show in their Fig. 21 radial profiles of the stellar velocity and velocity dispersion. The radial 
profile of the stellar velocity curve is very steep showing an increase in the stellar velocity 
from zero to 40-50 kms" 1 at 073 of the nucleus, which is in good agreement with our velocity 
curve along P.A.= 80° (although they do not show the two-dimensional velocity field). 



Davies et al. 



(120071 ) also report a decrease of a* from 130 kms 1 at 1"— 2" to only 70 kms 1 



at the very center. Our 0+ measurements shown in Figs.H] and [5] support the presence of a 



decrease in the velocity dispersion within the inner 2". But our values of 
down to ~ 125 kms -1 at 075 are larger than those reported by 



Davies et al. 



65 kms^ 1 at 2" 



(j2007l). Again 



we attribute this difference to different method s to obtain t 
well as to the distinct stellar templates used by IDavies et ai- 



r e velo city dispersion values as 



(120071 ) when compared to ours. 



Davies et al. 



(120071 ) have interpreted their observations of the rotation pattern plus low 



er* within the inner « 1" as being due to a compact stellar disk with « 70 pc radius that has 
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maintained the "cold" kinematics of the gas from which these stars have formed. They also 

claim that the stellar continuum increases by a factor of 2 above the inward extrapolation 

of the galaxy profile and that these two signatures are evidences of the presence of a nuclear 

disk which has experienced recent star formation. They have derived an inclination of 40° 

and a position angle of P.A.= 85° for the compact disk. This position angle is close to the 

major axis of the galaxy and the orientation and extent of the disk are consistent with the 

orientation and extent of the distribution of low a* values shown in our Fig.Hl 

l 



Another comparison we can make is with the results reported by lEmsellem et al 



( 120061 ). who have mapped the stellar velocity field in the inner ^20"of NGC 1068 using 
optical observations with SAURON, although with lower spatial resolution (apertures of 
f« 3"). They report a major axis position angle of 80°, derived from the stellar kinematics 
and which is also in agreement with the photometric major axis of the outer galaxy disk. 
This PA is consistent with our kinematics shown in Fig.Hl Their stellar velocity dispersion 
values are also similar to ours and they claim the presence of a a— drop within the inner 3", 
from values of ^150 km s -1 down to ~ 100 km s -1 . This is also in approximate agreement 
with our findings, only that we resolve better the "cr— drop", which we find to be within the 
inner arcsecond (70 pc). 



5.2. Stellar population 

5.2.1. Young stars and correlation with H2 ring 

Fig. [7] shows that the main stellar population components contributing to the flux 
within the inner 180 pc of NGC 1068 are the young and the young-intermediate ones. 
Within these components, we can verify from the output of the synthesis that the age 
bins which contribute with most of the flux are those of 30 Myr and 300 Myr, respectively. 



-24- 



Thus, it can be concluded that two main episodes of recent star formation have occurred 
in the nuclear region. The first burst occurred about 300 Myr ago, reaching as close as 0'.'5 
(35 pc) from the nucleus to the south-west, and north-east and north-west of the nucleus. 
The location of this burst may extend to even closer to the nucleus, but the dominance 
of the BB component in this region may preclude the detection of this component there. 
Particularly conspicuous is a region between 0'.'8 (60 pc) and 1'.'6 (100 pc) from the nucleus 
to the south-west, where it is the only stellar population component. The second burst 
occurred about 30 Myr ago, mainly in an approximately ring-like structure with an average 
radius of ~ 100 pc. 



In 



Vale et al. 



(120121 ) we have obtained emission line flux distributions from the same 



data we are using here. The flux in the molecular hydrogen emission line A2 . 1 2 1 8/xm, 



shown in the left panel of Fig. is distributed in a ring - first shown by 



(120091 ) - similarly to what is observed in NGC4151 (jStorchi-Bergmann et al 



Miiller- Sanchez et al. 



2010|). In 



Fig. [9], the H2 ring (left panel) can be compared with the young stellar population 
distributions x y and m y (central and right panels, respectively). This figure shows that 



and m y are enhanced at locations covered by the H2 ring, and, in particular, the region 
of strongest H2 emission is the region with the largest contribution of the young stellar 
component. 



The H? ring can be identified with the circumnuclear disk of molecular gas first 



reported by 



Schinnerer et al. 



20001 ). In Barbosa et al. (in preparation), we show that the 



H 2 emitting gas in this ring seems to be expandi ng; in the plane of the galaxy, confirming 



previous suggestions by 



Galliano & Alloin 


(2002 


) and results i 


the VLT ( 


Miiller-Sanchez et al. 


2009) 



20091 ). Recent observ ations of the 



molec ular gas emission in the mm spectral region (e.g. lines of CO, HCN) by iKrips et al. 



(120 111 ) also show this ring as well as its expansion, which, combined with the measured line 
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Fig. 9. — Comparison of the H 2 flux distribution (left) from 
tribution of the 30Myr stellar population percent contribution to the flux at 2.1 /im (central 
panel) and percent contribution by mass (right panel). The green contours on the central 
and right panels are from the H 2 flux distribution. 



H 2 A2.1218/xm 




arcsec 




Fig. 10. — Comparison of the continuum flux distribution (left panel) with the distribution 
of the percent contributi on of the BB com ponents (central panel) and the flux distribution 
of H 2 (right panel) from IVale et al.l ( 120121 ) . 
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ratios, have led these authors to propose a sce nario in whi c h par t of the molecular gas is 



radially blown outwards result of shocks. 



Krips et al 



(1201 if ) suggest that the shocks 



originate in the interaction of the AGN jet with the molecular gas. 



Krips et al. 



(120111 ) also point out that an expanding shell of shocked gas can originate 



in stellar winds from the evolution of a super stellar cluster 
favor this interpretation). This is the hypothesis put forth by 



although they do not 



Sakamoto et al 



(I2006h to 



explain the presence of "superbubble s" observed in molecula r gas emission in starburst 



galaxies such as NGC253 and M82 ( jMatsushita et al. 



20001 ). In the case of AGN, the 



presence of superbubbles has been rep orted by Lipari an d collaborators in a number of 



galaxies, example s being Mrk 231 



IRAS 04505-2958 (ILipari et al. 



e.g. 



Lipari et al. 



2005 



and references therein) and 



20091 ). These authors suggest that broad absorption line 



(BAL) systems seen in the spectra of AGN could be due to thes e superbubbles - g iant 



Lipari fc Terlevichl (120061 ) 



expanding gas shells which originate in circumnuclear starbursts. 
propose an evolutionary scenario in which massive circumnuclear starbursts in the vicinity 
of AGNs produce supernovae and hypernovae which lead to these large-scale expanding 
shells, which are often obscured by dust. Later in the evolution, as the starburst ages and 
the shells expand, the AGN becomes the main source of ionization of the circumnuclear 
gas. This scenario is shown to reproduce many of the observed correlations among 
spectral properties of in f rared- luminous AGNs. For lower luminosity Seyfert galaxies, 



Storchi-Bergmann et al. 



( 120011 ) also propose such an evolutionary scenario to explain the 
excess of young and intermediate age stellar population within the inner kiloparsec: the gas 
which flows in, first triggers star formation in the circumnuclear region, and the nuclear 
activity is then subsequently triggered. 

Thus, one possibility to be considered is that the H2 ring is actually a superbubble, 
for which the strongest emission comes from the galaxy plane, where there is more gas 
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to be swept by the bubble . This is sup p orted by the fact that NGC 1068 sits next to 



Mrk231 in the diagram of 



Davies et al 



( 120071 ) which relates the AGN luminosity with 



the circumnuclear starburst age. Although the expansion velocity of the H 2 ring - under 
100kms _1 (Barbosa et al., in preparation) - is much lower than in Mrk321 superbubbles, 
in nearby starbursts, observations in the milimetric spectral region show similar expansion 
velocities of a few tens of kms -1 . The morphology and sizes (100-200 pc) of superbubbles 
in nearby starbursts are also similar to those of the molecular ring of NGC 1068. One 
difference is the mass: while in the nearby starbursts the gas masses are of the order of 10 6 
solar masses (e.g. in NGC 253), the molecular g as mass in the inner palOOp c of NGC 1068 



is estimated to be an order of magnitude larger (jMuller-Sanchez et al. 
therein) . 



2009. and references 



Sakamoto et al. 



(120061 ) have 



In the case of the nearby starburst galaxy NGC 253, 
calculated the necessary power to drive the outflow and concluded that a 10 6 M super star 
cluster could be the origin of such bubbles, created by successive supernova explosions, 
which can have a duration of ~ 30 Myr, arguing that these bubbles are the end products of 
star formation. Having a mass an order of magnitude larger than the super-bubbles, it is 
not clear that the H2 ring in NGC 1068 has a similar origin. 

It is interesting to note that the apparent center of the H 2 ring does not coincide with 
the nucleus, but is located in the region w 1" to the south-west where the 300 Myr stellar 
population is particularly enhanced. If the H2 ring were actually a superbubble originated 
from stellar winds and supernova explosions, a cluster located there - in the center of the 
ring - could have originated it. A problem with this scenario is the fact that the ring 
is expanding and, even at a lower limit of the expansion velocity of lOkms -1 , it would 
reach its present radius of about 100 pc in only 10 Myr. As supernovae appear no later 
than pa 50 Myr, a cluster of age 300 Myr seems too old to have originated the H 2 ring in 
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NGC1068. 



And why would x y be enhanced in the H 2 ring? Why are the regions with strongest H 2 
emission (e.g. at ~ 1'.'5 east of the nucleus) the ones with highest concentration of young 
stars? Could an expanding bubble have triggered the 30 Myr episode of star formation? 
This seems unlikely due to the fact pointed out above that in 10 Myr the ring moves at 
least ~ 100 pc, and this is not consistent with the present radius of the ring, of ~ 100 pc and 
the age of the stars in the ring, of 30 Myr. 

On the other hand, the 30 Myr age is consistent with the presence of supernova 
explosions originated in this s ame population, which could be the origin of the shocks 



inferred by 



Krips et al. 



( 1201 ll ) from the molecular line ratios. The presence of such 
supernova explosions could also explain the high velocity dispersion observed in the region 
of strongest H 2 emission at 1'.'5 east of the nucleus (Barbosa et al. in preparation), and 
could be the origin of the apparent expansion of the H 2 ring. A spatial coi ncidence betwe en 
a 10-5 Myr stellar population and strong H 2 emission was also found by 



Martins et al. 



( I2010l ) at 100 pc south of the nucleus in a recent long-slit study of the stellar population of 



NGC1068. 



5.2.2. Older age components and continua 



A comparison between the stellar velocity dispersion (a*) distribution (FigHJ) and the 
distribution of x y i and Xj D in Fig. [7] reveals that, in most of the region where there is a 
drop in a*, the dominant stellar population is the intermediate age one (yi and io). We 



tentatively conclude that this intermedi ate age stel 



as we have found already for Mrk 1066 (IRiffel et al- 



ar po pulation is relat ed to the cr+ 



20101 ) and Mrk 1157 (IRiffel et al. 



Similar rings seem to be common in the central region of active galaxies (e.g. 



drop , 



20 



Barbosa 
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20061 ; iRiffel et al.l I2008L l2009al ) and can be interpreted as being due to colder regions with 
more recent star formation than the underlying bulge. This i nterpr etation is in agreement 



with the previous analysis of SINFONI data by 



Daviesetal 



(120071 ). who have estimated as 



200-3 00 Myr the age of the stars in the nuclear compact disk of NGC 1068. In a previous 



work, 



Thatte et al. 



(]1997af ) have claimed to have found a stellar core with intrinsic size 
scale of 45 pc. Using the estimated mass-to-light ratio and assuming some contamination 
from old stars, they propose an age of 500 Myr for a young age component in this region, in 
agreement with our results. 

We have found a small contribution from the old component, x , to the flux in the 
inner «0 / /5 (35 pc); but even though the flux contribution in the near-IR is not larger than 
~25%, the mass contribution is close to ~ 100%. The presence of an old stellar component 



aroun d the nucleus, with age > 10 9 yr, has already been suggested by 



Crenshaw fc Kraemer 



( I2000l ). on the basis of long-slit STIS spectra, as the dominant light component in the 



optical part of the spectrum. 

Fig.|H] shows the distribution of the power-law and black-body components. The 
power-law component does not contribute much at the nucleus, being spread around it. 
The vertical line is an artifact, which appears also in the lowest temperature black-body 
component; it is an instrumental effect, due to scattered nuclear light along the vertical 
mirror slices of NIFS. We attribute the extra- nuclear contribution of the power-law 
component to scattered light from the nucleus, which has been found to be a strong 



component of the continuum in the ultraviolet (I Crenshaw fe Kraemer 



2000|). 



As pointed out above, the black-body components dominate completely the light at the 
nucleus, in particular the ones with temperatures in the range 700-800 K (Fig. [8]), comprising 
almost 100% of the flux. We interpret them as ori ginating in the dusty torus surrounding 



the AGN, which is unresolved by our observations (IMarco fc AUoin 



2000 



Jaffe et al 



20041 ; 
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Raban et al. 



2009 



Martins et al 



2010(). In Fig.[10l we compare the continuum emission 



in the K band (left) with the sum of all black-body components (center) and with the H 2 
flux distribution (right). This comparison - which shows that the continuum emission and 
BB flux distributions are similarly concentrated - gives support to the identification of 
the peak of the K-band continuum with the black-body components contribution. Some 
contribution of the black-body components is seen also beyond the inner arcsecond, which 
we attribute to the presence of faint dust emission or to further contribution from scattered 
nuclear light spread in the nuclear region. The comparison between the central and right 
panels of Fig. [TU] shows that some BB emission is seen also along the west part of the ring 
(to the right in Fig. [TO]) . 



6. Conclusions 

We have mapped the stellar kinematics and age components, as well as the contribution 
of power-law and black-body emission components in the inner 180 pc radius of NGC 1068 
at a spatial resolution of 8pc using integral field spectroscopy in the near-infrared. The 
main results of this paper are: 

• The stellar kinematics are well described by rotation with a projected amplitude of 
i^45kms _1 and major axis a position angle of ~80deg; 

• The stellar velocity dispersion a* shows a drop from 165 km s" 1 to 125kms _1 within 
the inner 1"- 70 pc at the galaxy; 

• The dominant stellar populations over the inner 180 pc radius are the young (mainly 
30 Myr) and intermediate age (mainly 300 Myr) ones; 

• Within the inner 0'.'5 (35 pc), the black-body components with 700<T<800K 
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dominate the flux, and are attributed to the dusty torus (unresolved by our 
observations) surrounding the active nucleus; 

• Within the inner 35 pc, the old component (age> 5 Gyr) dominates the contribution 
of the stellar population in terms of mass fraction; 

• The distribution of the youngest (30Myr) stellar population is co-spatial with an 
expanding ring of H 2 emission at ^lOOpc from the nucleus. We speculate that 
supernova explosions from this young stellar population could be the origin of the 
apparent expansion of the ring; 

• In the region of the cr*-drop, the dominant stellar population is the intermediate age 
(mainly 300 Myr) one. 



NGC 1068 - the best studied Seyfert 2 galaxy is thus another case of an active nucleus 
surrounded by rece nt star formation in a circumnuclear ring, supporting the evolutionar y 



scena rio for AGN (IStorchi-Bergmann et al. 



2001 



Lipari fc Terlevich 



2006 



Davies et al. 



20071 ). In such a scenario, gas from the external regions of the galaxy (or even acquired 



from outside) flows to the nuclear region, first triggering star formation in the inner few 
hundred parsecs, and subsequently reaches the nucleus to feed the supermassive black hole, 
triggering the nuclear activity. 
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